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In a recent paper we have considered an extended Skyrme model which includes fourth and 
sixth-order terms and have investigated the soliton breathing mode in the framework of the linear 
response theory. We have found that the monopole response function exhibits an unbound sharp peak 
which we have identified to the Roper resonance A^(1440). Since the fourth and sixth-order terms can 
be derived from a local approximation of an effective model with p Q and uj mesons Q respectively, 
it seems therefore interesting to explore the soliton breathing mode, with the same method as in 
[Ql (first proposed in Ref. ||5|), within an effective Lagrangian which incorporates these two mesons 
explicitely. This is the purpose of the present paper. 

1. Our starting point is the following Lagrangian density: 

C = IItt {d^Ud^U^) + C^p + C^^ + -jlmlTiiU + -2) (1) 
where U is an SU(2) matrix parametrized by the pion fields tTq, normalized to the pion decay constant 

= exp[ m.T/U ] , (2) 

the Tq's being the usual Pauli matrices. The first term in Eq. (||) corresponds to the well known 
nonlinear a model. The second one corresponds to the coupling between the isospin one p-meson and 
the pion field. The p vector meson is frenquently introduced in the chiral Lagrangian as hidden gauge 
particle or as a massive Yang-Mills field In both cases the soliton was shown to be unstable 
and consequently one can not describe the excited states, in particular the breathing mode. In a 
recent paper it was proven that the soliton is stable with respect to the breathing fluctuations if 
the p-meson is introduced in the chiral Lagrangian following the standard prescription of nonlinear 
chiral symmetry for massive particles In Ref. however, the p-meson is described in terms of 
antisymmetric tensor field as proposed by the authors of Refs. |10, 11]. It has been shown in Ref. [12 



that this description is canonically equivalent to a vector field formulation provided the Skyrme term 
[13| is added to the latter in order to keep the corresponding Hamiltonian bounded. In this work we 



make use of this last formulation to introduce the p-meson. It yields: 

C^^ = ^Tr{V,Vn-lTrl^{V^V,-V,V^ + i^[up,u,] )^ | (3) 
where is the p-meson field and 

= + [r^, ], = ^[u\d^u], u = = m^d^^Uu^ . (4) 

The two constants appearing in Eq. (^) are the p-meson mass Mp and the dimensionless parameter 
gv which can be related to the p — > vrvr decay width. 
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The third term in Eq. corresponds to w-meson exchange It reads: 

J^TTio = -\i ^tii^u - d^ujf, f + ^mluj^uj'' + P^ujf.B'' (5) 

where is the omega field and i?^ the baryon current |p!3| , |l4| 

= ^e'^'^^Tr { {d,U)U\daU)U\df,U)U^ }. (6) 

The two new constants appearing in Eq. (^) are the w-meson mass m^^ and the parameter which 
can be related to the uj — > tt'j width |15]. Finally, the last term in the Lagrangian (jl]) implements a 



small explicit breaking of chiral symmetry, m-j^ being the pion mass. 

2 . In order to work with true degrees of freedom and eliminate the constrained fields Vq and 
u!q (the corresponding canonical momenta conjugate are equal to zero), it is suitable to consider the 
Hamiltonian density. Let and being the canonical momenta conjugate to the pion field vf, the 

p field and the uj field respectively. We make the hedgehog ansatz for the pion and the most general 
spherical ansatz for the p and the uj fields: 

7f = /^Fr, Vi = U{vi{Ti - {T.r)fi) + V2{T.r)fi - V3{t X r)i} , uJi = ^uj fi, 

- /2 (7) 
$ = -JL(j) f , Pi = — {71 (Ti - (r.f)ri) + 72(T.f)fi - 73(f x r)i} , pi = —p fi . 

The energy functional expressed in terms of the time-dependent radial functions F, (p, Vi,^i,u:,p reads: 

fOO 111 1 

E = 47r/2 dr{ -7? + -7I + -7I + My{2vl + v^,) + -mlrW+ 



l{ct> + 25-73 - P-ujf + 2{{rviy - V2cf + Tv^(27ic - (720')'+ 
z r r AMpT'^ 

1,2j7/2 , „2 I r,),^2„2„,2 I 1„2 , _2„2 



(8) 



-r^F'^ + + 2Myv^ + -p^ + mir\l - c) + 
^2 ^1 

{2v3C + g-f + 2igsF' + (rvsYf + .^^{Ps'F' + {rpYf } 
r zm^r 

with g = \/2gv/ fn, jS = f^uj/'^'^'^fiT, s = sinF and c = cosF. Primes indicate radial differentiations. 
By solving the static Hamilton equations, i.e., 

F = 4> = Vi='yi=uj=p = 

where dots indicate partial time differenciations, we find that the only nonsingular solution which 
extremizes E is the one which corresponds to 

(j)(r) = = = 71 = 72 = 73 = u; = . (9) 
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The set of the three other static Hamilton equations corresponding to F, vs and p reads: 

s — 
r{rF)" + 4g--{Mpr^ + s'^)v3 + 2sc{4vl - 1) + l3—p - m^r^s = , 

r{rv3)" - (2c2 + My)v3 - g(,— + rc") = , (10) 

_ .2 

rV - (2 + mlr'^)p + r^(3{-^F'y = . 

We solve numerically the equations ( [l0| ) with the boundary conditions F(0) = vr, F(oo) = in order 
to find a winding number one solution. In Fig. 1 we plot that solution for the following set of physical 
parameters (fixed by fitting to the low energy meson observables) : 

U = 93, Mp = 769, = 782, m„ = 139.5 (MeV) 
= 9.3, gv = 0.09 . 

The classical soliton mass we find is M = 1515 MeV, and is lower than the mass which corresponds 
to the local approximation (1714 MeV) that one finds if the masses Mp and and the coupling 
constant Pui are increased to infinity, keeping gy and the ratio Puj/muj finites. We also plot in Fig. 1 
the local approximation's chiral function for reference. 

3. In order to investigate the soliton breathing mode within the model (||) we consider the 
linearized time-dependent Hamilton equations in presence of an external infinitesimal monopole field 
with a frequency Q. This can be done by adding the following perturbation term [Q, ^ 

<5i?int = Mm exp(r/t) (12) 

to the energy functional (P). In Eq. ([l2|), (r^) is the isoscalar mean square radius [^] and r/ a 
vanishingly small positive number. 

Because the term ( p^ ) is weak it introduces only small changes of the classical Hamilton solution. 
We thus can use the linear approximation and look for a solution in the form 

F{r, t) = F{r) + e[ 6F{r)e'^^-''''^^ - c.c ]/2i , 0(r, t) = 0(r) + e[ 6^{r)e'^^-'''^^ - c.c ]/2i , 
oj{r, t) = u{r) + e[ 5oj{r)e'^^-''^'^^ - c.c ]/2i , p{r, t) = p{r) + e[ 5p(r)e*(^-^'')* - c.c ]/2i , 
Vi{r,t)=v^{r) + e[5vi{r)e'^''-'^^' -c.c]/2i , 7,(r, t) = 7,(r) + e[ 57i(r)e*(^-^'')* - c.c ]/2i , 

with the boundary conditions that all the fiuctuations and their time derivatives are equal to zero at 
t = —oo. We are not going to write all the ten linearized equations but it is straightforward to show 
from Eq. (^) that the equations corresponding to the fluctuations 5(j),5uj,5^i^5vi and 5v2 decouple 
from those which correspond to 8F^ Sv^, and 8p, and consequently do not contribute to the breathing 



4 



(13) 



6^ 



.0, 



mode. The equations corresponding to the breathing fluctuations 5F, 5v3 and 6p read : 

/All Ai2 

A21 A22 
\^3i A33, 

where Sip and 5^ are combinations of 6vs and 6p, namely 

Stp = 6v3 + g-6F , 5C = 6p + j3—5F . 

r r 

The matrix elements Aij of the operator A appearing in Eq. ( p!^ ) are themselves operators and read 
as: ^ 

^11 = —r + — [ {2vl - -){s^ - c^) + -s-{-s~ + V3C) + -ivi;r~(-s~ - V3C) 



1 



/5 



^12 = -J[^(^' + M,V)+4^3C] 



s(2cp s ) H — m r c ] 

4r r 4 



A 
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«^ 5. 2 



^2i = -^[-(s' + M2r2)+4t;3C 



(14) 



^22 
^31 
^33 



-mlf- 



We recall that s means sin(i^), c stands for cos(i^) and F is the static solution for the pion field. 

The monopole response function is determined from the evolution of the isoscalar mean square 
radius of the soliton with respect to the frequency fi. The mean square radius is given by [16| 

dV5°(r,t)r2 = -_/ dr r2sin2(F)F' 
Jo 

where B^{r,t) is the time-component of the baryon current (P). Up to first order in e, (r^) reads 



(r2)(i) = (r2)o + i:[/(0)e^M* 



2i' 



c.c 



where / is the linear response function 
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f{n) = - dr rsm^{F)5F{r). 
vr Jo 



(15) 



The spectral representation of the response function / can be extracted by using equations (13). It 
reads 



vr 



(16) 



where the hmit ij ^ is, as usual, imphcit, and corresponds to the boundary conditions specified 
above. In Eq.(|l6|) the state x is defined by 

{X\r) = ^-(1 ,0,0) (17) 



and the Xn are the eigenstates of the operator A (see Eqs. (13), ([T^), with the eigenvalues A^, 
normalized according to 

iXnlXm) = / /|?-^dr Xnir) Xmir) = 5nm ■ 

Jo 

Since the distribution of collective strength is directly related to the imaginary part of the linear 
response function, we should consider this quantity. We display in Fig. 2 the imaginary part of the 
function (|l6| ) with respect to the frequency Q for the parameters given in Eq. (pA[). We find that it 
exhibits a pronounced peak at 400 MeV which we identify with the excitation energy of the Roper 
resonance. In Ref. ||l| we have found nearly the same value within the local approximation to the 
Lagrangian (|^) except for the detail that the dimensionless Skyrme term parameter e [jl3[ taken in ||l| 
is equal to 7.1 [0, [l^, while the value of e which corresponds to the parameters (^) we use in this 
work is e = l/2gv = 5.57. We show in Fig. 2 the response function which corresponds to the local 
approximation taking e = 5.57. One sees from this figure that there is an improvement over the 
local case, although a slight one (~ 40 MeV) but in the good direction. A more realistic improvement 
consists of extending the model (|l]) to include other mesons such as the vector meson-^i and the 
scalar meson-e despite of the increasing complexity implied. 

4. To summarize, we have investigated the soliton breathing mode within an effective La- 
grangian with finite mass mesons (the p- and the w-mesons) in the framework of linear response 
theory. Our results for the location of the Roper resonance show an improvement over the local 
approximation where these two mesons are taken to be infinitely heavy. 

It is quite well established (see for example [|15|) that for the description of static nucleon properties, 
finite mass meson Lagrangians are more accurate that their local residues, such as the Skyrme model. 
As this work suggests, this feature is also true for the description of dynamical observables such as 
the breathing mode. 

Finally, the difficulty to find the Pll resonance A^(1440) in analysis based on phase shifts method 
[19 1 and the encouraging results obtained in this work and in Refs. ||l], ^ lead us to think that the 



linear response theory is the most adapted method to describe the Roper resonance within effective 
Lagrangians. 

We are grateful to J. Aichelin, D. Kalafatis, B. Moussallam and D. Vautherin for a critical reading 
of the manuscript. 
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Figure captions 



FIG. 1. The solution of Eq. ( p!0| ) minimizing the energy functional (^). The full line displays the 
static function F{r) in this model, while the dotted one corresponds to the local approximation of (||). 
In the dashed and dashed-dotted lines we display the p and the uj degrees of freedom respectively. 
Namely f3(r) and p{r). 

FIG. 2. Imaginary part of the response function / (fm^) versus the energy $7 in MeV. The full 
line corresponds to the Lagrangian with the parameters given in Eq. (11) while the dotted line 
corresponds to the local approximation of with f^^ = 93 MeV, fS^j = 9.3 and e = l/2gv = 5.57 (c.f. 
Ref. M). 
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